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Changes in zonal surface temperature gradients and Walker circulations in a wide range 

of climates 
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ABSTRACT 

Variations in zonal surface temperature gradients and zonally asymmetric tropical overturning circu- 
lations (Walker circulations) are examined over a wide range of climates simulated with an idealized 
atmospheric general circulation model (GCM). The asymmetry in the tropical climate is generated by 
an imposed ocean energy flux, which does not vary with climate. The range of climates is simulated 
by modifying the optical thickness of an idealized longwave absorber (representing greenhouse gases). 

The zonal surface temperature gradient in low latitudes generally decreases as the climate warms in 
the idealized GCM simulations. A scaling relationship based on a two-term balance in the surface 
energy budget accounts for the changes in the zonally asymmetric component of the GCM-simulated 
surface temperature gradients. 

The Walker circulation weakens as the climate warms in the idealized simulations, as it does in 
comprehensive simulations of climate change. The wide range of climates allows a systematic 
test of energetic arguments that have been proposed to account for these changes in the tropical 
circulation. The analysis shows that a scaling estimate based on changes in the hydrological 
cycle (precipitation rate and saturation specific humidity) accounts for the simulated changes 
in the Walker circulation. However, it must be evaluated locally, with local precipitation 
rates. If global-mean quantities are used, the scaling estimate does not generally account for 
changes in the Walker circulation, and the extent to which it does is the result of compensat- 
ing errors in changes in precipitation and saturation specific humidity that enter the scaling estimate. 



1. Introduction 

The evolution of zonal surface temperature gradients 
in the tropics is of interest in the study of future and 
past climates clima tes (e.g., Knutson and Manabd 119951 
iFedorov et al. 2006 ). Recently, the changes in the zonally 
asymmetric component of the tropical overturning circu- 
lation, the Walker circulation, have been highlighted as 
a robust response to warming in climate change s i mula- 
tions with comprehens ive GCMs ( Held and Sodenl 2006 ; 
IVecchi and Soden 2007 ). Here, we examine the response of 
these zonally asymmetric aspects of the mean tropical cli- 
mate to radiative perturbations in an idealized GCM and 
test theoretical arguments designed to account for their 
changes. 

How zonal surface temperature gradients evolve as the 
climate warms is unclear. On one hand, increased up- 
welling of cold water on the eastern ocean margins may 
lead to inc reased zonal tempera ture gradients in transient 
warmings (jClement et all 1 19961 ). and this mechan ism ap- 
pears to be important in coupled GCM simulations (jHeld et al 
20101) . On the other hand, paleoclimatc evidence from the 



Pliocene suggests substantially weakened or collapsed zonal 
temperature gradients in a climate o nly somewhat warme r 



than Earth's current c l imate (e.g., IFedorov et al.l 120061 ) . 
(However, see Wunsch ( 20091 ) for an alternative interpre- 
tation of the data.) The multi-model mean of the fourth 
Intergovernmental Panel on Climate Change (IPCC) as- 
sessment models shows little change in zonal temperature 
gradients in the Pac ific basin over the twenty-first century 
(DiN ezio et al ]l2009l) . As a step towards understanding the 
full climate dynamics that control the zonal surface tem- 
perature gradient in the tropics, we develop a scaling rela- 
tionship based on the surface energy budget and compare 
it to the results of the idealized GCM simulations. 

The connection between changes in the hydrological 
cycle and large-scale atmospheric circulations in response 
to radiative forcing, for instance, by changed greenhouse 
gas concentrations, can be subtle. For example, the mass 
fluxes in the Hadley cell and extratropical eddy kinetic en- 
ergy can vary non-monotonically as a longwave absorber is 
systematically varied in an idealized GCM, with implica- 
tions for water vapor flu xes and the distribution of evapor a 



tion and precipitation ( O'Gorman and Schneider 2008b. a 
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Schneider et al. 2010t Levine and Schneider 2011 ). 

Analysis of comprehensive GCM simulations suggests 
the connection between changes in the h ydrological cycle 
and upward mass fluxes is stra ightforward ( Held and Sodenl 
20061 Ivecchi and Sodenll2007h . The argument is that pre- 



cipitation, which is constrained by the energy balance at 
the surface in the global mean, increases with warming 
more slowly than near-surface specific humidity, which in- 
creases similarly to the rate given by the Clausius-Clapeyron 
relation for small changes in relative humidity; thus, the 
upward vertical velocity that gives rise to condensation of 
water vapor weakens. This scaling argument is sometimes 
interpreted to account for cha nges in the Walker circu- 
lation (e.g., IVecchi et al. 2006 ). in part because changes 
in the zonally asymmetric component of the circulation 
are larger than those i n the zonally symmetric component 
( Held and Sodenll2006l) . A difficulty in applying the hydro- 
logical cycle argument directly to local circulations, such as 
the Walker circulation, is that the surface energy balance 
constraint on the precipitation is lost, as water vapor con- 
vergence depends locally on the atmospheric circulation. A 
central goal of this work is to assess the extent to which the 
hydrological cycle arguments apply locally to the Walker 
circulation, and to compare this to global-mean versions 
of the scaling arguments to determine the extent to which 
they are adequate. 

Here, we use an idealized GCM to systematically exam- 
ine the response of the zonally asymmetric tropical surface 
temperature and the Walker circulation to changes in the 
atmosphere longwave opacity. A scaling estimate based 
on the surface energy balance accounts for the changes in 
the zonal surface temperature gradient. We assess the hy- 
drological cycle estimate for changes in the Walker circu- 
lation by applying it locally. The scaling estimate is ac- 
curate, though it requires knowledge of the precipitation 
field. In contrast, the global-mean version of the scaling 
estimate does not in general account for the Walker circu- 
lation changes. 

2. Idealized GCM 

a. Model description 



scribed in O'Gorman and Schneider (l2008bh . There is no 
condensed water phase in the atmosphere (i.e., there are 
no clouds). The climate is perturbed by varying the opti- 
cal depth in the gray radiation scheme by a multiplicative 
factor to mimic variations in the trapping of longwave ra- 
diation due to variations in greenhouse gas concentrations. 
The primitive equations are integrated using the spectral 
transform method in the horizontal with T85 resolution 
(but the simulation results arc similar with T42 resolution) . 
Finite differences in a — p/p s (with pressure p and surface 
pressure p s ) discretize the vertical coordinate, with 30 un- 
evenly spaced levels. The time-stepping is semi-implicit, 
with a leap frog timestep of 300 s. The simulation statistics 
are averaged over the last 1100 days of the simulation fol- 
lowing a spinup period of either 700 days for an isothermal, 
resting initial condition or 300 days for an initial condition 
from an equilibrated simulation with similar optical thick- 
ness. Aside from the horizontal resolution and timestep, 
the only difference betw een the simulations presen t ed here 
and those presented in lO'Gorman and Schneider! (|2008bl ) 
is that we include a representation of ocean energy flux 
divergence in the slab ocean and present fewer simulated 
climates. 

b. Prescribed ocean energy flux divergence 

The surface boundary condition is a slab ocean that ac- 
counts for the radiative and turbulent surface fluxes. Ad- 
ditionally, it has a prescribed ocean energy flux divergence 
that represents the influence of the ocean circulation on 
the surface temperature. This is commonly referred to as 
a Q-flux in climate modelling. 

The Q-flux consists of a component that varies only in 
latitude and represents the divergence of the meridional 
ocean energy flux 



where <j> is latitude, 



4>l 



CXP 1 "ft 



(i) 



16°, and Q n = 50Wm- 2 fo l 



We use an idealized, moist primitiv e-equation atmo- 
spheric GCM similar to those descr ibed in lO'Gorman and Schn 
(|2008bl) and iFrierson et al.l (j2006l) . Briefly, the model in- 



lowing Bordoni ( 2007 ) and Bordoni and Schneider ( 20081) . 
The meridional component of the Q-flux is necessary to 
have an Earth-like tropical circulation in the reference cli- 
mate; without it, the zonal-mean (Hadley) circulation is 
icntly strong that there is mean ascending motion ev- 



cludes an active hydrologi cal cycle with a simplified Betts- 
Miller convection scheme (|Friersonll2007h . The scheme re- 
laxes over a fixed timescale convectively unstable columns 
to moist pseudoadiabatic temperatures with constant (70%) 
relative humidity. It has a gray longwave radiation scheme, 
and a slab-ocean surface boundary condition with the heat 
capacity of 1 m of water. The top-of-atmosphcre insolation 
is an idealized, annual-mean profile, and there is shortwave 
absorption in the atmosphere, both of which are as de- 



ery where along the equator (i.e., the descending branch of 
the Walker circulation is weaker than the Hadley circula- 
tion). 

Zonal symmetry is broken by adding and subtracting 
Gaussian lobes along the equator that are spaced 180° lon- 
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Fig. 1. Ocean energy flux divergence, Q-flux, with positive values indicate cooling tendencies. Large panel: longitude- 
latitude Q-flux contours with contour interval of 20 Wm -2 and negative values are dashed. The gray rectangles indicate 
the averaging regions. Right panel: latitude vs. zonal-mean Q-flux. Bottom panel: equatorial Q-flux (averaged within 8° 
of the equator) vs. longitude. 
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(0.6, 0.7, 0.8, 0.9, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.5, 3.0, 4.0, 
6.0). The range of optical depths presented here is re- 
duced compared to lO'Gorman and Schneiderl (|2008bf ) : sim- 
ulations with the smallest optical depths are omitted be- 
cause the fixed meridional Q-flux leads to reversed Hadley 
cells and concomitant double intertropical convergence zones 
in sufficiently cold climates. In total, fourteen simulations 
are presented here. 



where A is longitude, Ai = 30°, Xe 
0! = 7°, and Qj = 40 Wm" 2 . 

The resulting Q-flux, shown in Fig. [TJ has simil ar mag- 



270°, 3. Zonal surface temperature gradients 



nitude and spatial scale to reanalysis estimates (e.g.. lTrenberth 
200 ll their Fig. 5). 



a. Simulation results 
et af.l gure[3]shows the simulated zonal surface temperature 



c. Series of simulations 

As in O'Gorman and Schneider ( 2008bl ). the longwave 
optical depth is varied to mimic changes in greenhouse gas 
concentrations. The optical depth r = ar rc f is the product 
of a multiplicative factor a and a reference profile r re f that 
depends on latitude and pressure, 

T re f = [ficr + (1 - fi)cr A ][T eq + (r p - T eg )sin 2 </>], (3) 



where fi = 0.2 and the reference optical depths at the 
equator and pole arc r eq = 7.2 and t p = 1.8, respec- 
tively. We conduct simulations with rcscaling factors a = 



differences averaged over the regions indicated in Fig. [T] for 
each of the equilibrated climates, which arc identified by 
the time- and global- mean (denoted (•)) surface temper- 
ature (T s ). Following 15'Gorman and Schneiderl (|2008bl ). 
the reference climate with a = 1.0 is indicated by a filled 
symbol. The east-west surface temperature difference AT, 
decreases nearly monotonically as the climate warms and 
varies by more than a factor of three over the simulated 
range of climates. 

The simulated zonal temperature difference in the ref- 
erence climate (1.4 K) is smaller than that in Earth's cur- 
rent climate (~5 K). While we have used a Q-flux that has 
comparable magnitude to estimates for the current climate, 
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Fig. 2. Zonal surface temperature difference (west mi- 
nus east in Fig. [1]) vs. global-mean surface temperature 
of GCM simulations (circles with black dashed line), and 
scaling estimate (|8| (gray line) . The averaging conventions 
used to evaluate the scaling estimate are described in the 
text. 



the lack of clouds and the concomitant zonal asymmetry 
in cloud radiative forcing leads to a smaller zonal asymme- 
try in the equilibrated surface temperature. One remedy 
would be to impose a larger asymmetry in the Q-flux, since 
in quasi-cquilibrium theory, there is a near equivalence 



betw een surface fluxes and radiative forcing (jSobel et al 



201Ct section 2.3). However, we chose not to amplify the 
Q-flux asymmetry. Therefore, it is more meaningful to 
consider fractional changes in the simulated zonal surface 
temperature contrast than the magnitude of the simulated 
changes. Our results hence may give changes in the sur- 
face temperature contrasts that, in reality, would likely be 
modified by processes we ignored, such as cloud-radiative 
feedbacks. 

b. Scaling estimate 

Figure [3] shows the surface energy budget as a function 
of longitude for the Earth-like reference simulation. Ex- 
amining the variations in longitude of the surface energy 
balance reveals that the dominant balance in the zonally 
asymmetric component is between the prescribed Q-flux 
and the evaporative fluxes. This is to be expected in Earth- 
like and warmer climates where the Bowen ratio (ratio of 
sensible to latent surface fluxes) is small in the tropics. An 
analogous ratio can be formed for the ne t longwave radia- 
tion at the surface ( Pierrehumbert 2010h . and this is also 
small for sufficiently warm climates. 

This two-term balance between the zonally asymmetric 
component of the Q-flux V ■ F\ and evaporation holds over 
the range of climates we simulate, so that 



where A refers to east-west differences. 

As the Q-flux (ocean energy flux divergence) is fixed, 
the bulk aerodynamic formula for evaporative fluxes can 
be manipulated with suitable assumptions to invert the 
zonal difference in the evaporative fluxes into a zonal dif- 
ference of the surface temperature. The evaporative fluxes 
are computed by the bulk aerodynamic formula 

E = pca||v||( 9a (T.) - q) « pc d ||v||(l - H)q s (T s ), (5) 

with air density p, drag coefficient Cd, surface wind speed 
||v||, saturation specific humidity q s , specific humidity of 
surface air q, surface temperature T s , and surface relative 
humidity TL. The second formula is approximate in that 
the ratio of the specific humidity and the saturation specific 
humidity is set equal to the relative humidity and in that 
the air-sea temperature difference is neglected. Neglecting 
the air-sea temperature difference is justifiable because the 
subsaturation of air dominates over the air-sea temperature 
difference for conditions typical of tropical oceans. 

For the climate changes we are examining, the east- 
west evaporation difference then scales like the east-west 
difference in the surface saturation specific humiditjQ 



AE ~Aq s (T s ). 



(6) 



This is the result of the strong dependence of the saturation 
specific humidity q s on temperature (fractional changes 
of ~7%K _1 for Earth-like temperatures). In contrast, 
changes in the surfa ce relative humidity H are relatively 
small. For example, ISchneider et al. (2010) presented an 
argument based on a simplified surface energy budget that 
suggests ~1% K _1 changes in H. Regional changes that af- 
fect the zonal gradient of relative humidity are potentially 
important, but neither comp rehensive GCM simulations 
(jCGorman and Mullen |201 Of) nor the idealized GCM sim- 



ulations presented here exhibit substantial changes in the 
zonal relative humidity gradient. The neglected changes in 
the air-sea temperature difference and surface wind speed 
are also generally smaller than the changes in the surface 
saturation specific humidity. 

To obtain an cast-west surface temperature difference, 
rather than a saturation specific humidity difference, we 
linearize the saturation specific humidity 



Aq s (T s 



dT 



-AT, 



(7) 



The resulting scaling estimate for the east-west tempera- 
ture difference is the product of the east-west saturation 
specific humidity difference Aq s and the inverse of the rate 



A (V • Fx) ~ ALE, 



(4) 



We are not suggesting that this scaling holds in general (see, 
e.g., O'Gorman and Schneider (2008b) and our Fig. [6] for a coun- 
terexample in the global mean), but only for the zonal difference in 
the tropics. 
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Fig. 3. Surface energy fluxes of the reference simulation vs. longitude averaged within 8° of equator, "lw" and "sw' 
refer to the net longwave and shortwave radiative fluxes at the surface, respectively. 



of change of the saturation specific humidity with temper- 
ature: 



(8) 



For invariant ocean energy fluxes, as in the idealized GCM 
simulations, Aq s is approximately constant by Eqns. (j4]) 
and and the east- west surface temperature difference 
will decrease with warming by an amount given by the rate 
of change of the saturation specific humidity with temper- 
ature (i.e., given by the Clausius-Clapeyron relation as the 
surface pressure is nearly constant), evaluated using a tem- 
perature that is representative of the surface in the tropics. 

The decrease in the zonal temperature difference in a 
warming scenario can be anticipated by considering the 
counterexample of invariant zonal surface temperature gra- 
dients: If the east-west temperature difference remained 
the same as the mean climate warmed, in the absence of 
changes in near-surface relative humidity, the east-west dif- 
ference in evaporative fluxes would rapidly increase with 
warming. The resulting zonal asymmetry in the evapo- 
ration would be too large to be balanced by an invariant 
ocean energy flux, and the other terms in the surface energy 
balance are too small to restore equilibrium. In a transient 
adjustment, the excess in evaporation would lead to a re- 
duction in the east-west surface temperature asymmetry 
as, for example, the large evaporation over warm regions 
would cool the surface. The sensitivity of the evapora- 



tive fluxes to climate change and concomitant reduction in 
zonal surface temperature gradients with warming was doc- 
ument ed in the GCM simulations of Knutson and Manabei 
(1995), who called it "evaporative damping" . 

These approximations are adequate for the longwave 
radiation-induced climate changes in the idealized GCM 
simulations, though they may not be in genera l. For e x- 



ample , in analyses of climate change simulations. Kie et al 



( 201Cf ) showed that changes in surface winds must be ac 
cou nted for in the surfac e energy budget of some regions, 
and DiNezio et al. ( 2009() showed that cloud-radiative feed- 
backs affect the net surface shortwave radiation. 

c. Assessment of scaling estimate 

Fig. [U shows the scaling estimate ([5]) in gray. The 
derivative of the saturation specific humidity with respect 
to temperature in (|SJ) is evaluated with the same simpli- 
fied formulation for the saturation specific humidity as is 
used in the GCM, using the simulated time- and zonal- 
mean surface temperature averaged within 8° of latitude 
of the equator. The saturation specific humidity differ- 
ence Aq s , evaluated over the averaging regions in Fig. [1] 
varies by ~15% over the range of simulations, and we use 
the mean value, 2.4 x 10 -3 , in the scaling estimate. The 
magnitude of the variations in Aq s is consistent with the 
magnitude of the neglected terms in the bulk aerodynamic 
formula. Overall, there is close agreement (maximum devi- 
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Fig. 4. Contours of the zonally asymmetric component of 
the time-mean pressure velocity uj* in the longitudc-sigma 
plane averaged within 8° of the equator for three simula- 
tions with global- mean surface temperature (T s ) indicated 
on the right. The contour interval is 0.01 Pas -1 with posi- 
tive values (descending motion) dashed. The zero contour 
is omitted. 



ations ~15%) between the scaling estimate and the results 
of the GCM simulations. 

We note that while the simulations and scaling estimate 
feature rapidly decreasing zonal temperature gradients as 
the mean temperature increases (7.4% K _1 at the reference 
climate), the decreases are not as rapid as proxies of the 
Pliocene suggest (~15-30% K _1 ). But it is possible that 
the neglected radiative and ocean feedbacks increase the 
sensitivity to changes. 



a. Simulation results 

Figure [4] shows the zonally asymmetric component of 
the time-mean pressure velocity, Dp/ Dt = uj, in the longitudc- 
sigma plane for three simulationtQ The zonally asymmet- 
ric pressure velocity is large over the regions of perturbed 
Q-flux; it is small elsewhere. The pressure depth over which 
it has large amplitude increases with warming because the 
tropopausc height increases; its magnitude decreases with 
warming. Note that a consequence of the Walker circu- 
lation's localization over the Q-flux perturbations is that 
the averaging convention we have chosen is adequate to 
capture nearly all of the ascending and descending mass 
fluxes — it is a closed circulation and alternative measures, 
such as surface pressure differences, yield similar results. 

Figure [5] shows that the Walker circulation mass flux in 
the ascending region rapidly weakens with warming. The 
fractional rate of decrease is 4.4% K" 1 at the reference 
climate relative to the global-mean surface temperature. 
(The fractional change is larger relative to local surface 
temperature changes as the warming is greater in high lat- 
itudes.) The Walker circulation varies by a factor of six 
over the simulated range of climates. In this figure, the 
zonally asymmetric component of the pressure velocity is 
evaluated on the sigma surface where it is maximum, which 
varies from a = 0.60 in the coldest climate to a = 0.32 in 
the warmest climate. These changes roughly follow those 
of the zonal-mean tropical tropopausc (averaged within 8° 
of the equator) which rises from a — 0.24 to a = 0.03 deter- 
mined by the World Meteorological Organization (WMO) 
2Kkm~ 1 lapse rate criterion. If the Walker circulation is 
evaluated using fixed sigma levels in the lower troposphere, 
the weakening with warming is more rapid owing to the 
combined changes in the vertical structure of the circula- 
tion and the changes, shown in Fig.[5j that are independent 
of the vertical structure. 

b. Scaling estimate 

In saturated ascending motion, the vertical advection 
of saturation specific humidity q s leads to condensation c, 



/a 



C, 



(9) 



where ui^ = u>H(—u>) is a truncated upwar d pressure veloc 



ity a nd H is the Heaviside step function ([Schneider et al 
2010h . For a mass- weighted vertical average taken from 



the lifting condensation level to the tropopause (denoted 
by {■}), the precipitation P balances the vertical advection 
of specific humidity: 



{w%q a }KiP. 



(10) 



2 A11 figures show p 3 Dcr/ Dt, which differs only slightly (~ 0.1% in 
the region of interest) from Dp/Dt as crDp s /Dt is small. 
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Fig. 5. Left: zonally asymmetric component of the pressure velocity to* (black) and zonally asymmetric component of 
the upward pressure velocity uj^* (red) averaged over the ascending region of the Walker circulation (cf. Fig. [1]). The 
pressure velocity oj* is evaluated on the sigma level where it is minimum (strongest ascent), and oj^* is evaluated on 
the same level. The upward velocity has been multiplied by a factor of 1.6 determined by minimizing the mean squared 
deviation between u^* and u*. Right: zonally asymmetric component of the pressure velocity, w*, (black) and hydrological 
cycle scaling estimates: —gP*/A v q s (blue), — gP*/[q s b] (magenta), and —g(P)/(q s b) (orange). The scaling estimates are 
multiplied by constants of 1.5, 1.7, and 0.8, respectively, chosen to minimize the mean squared deviations between the 
scaling estimates and a;*. The near-surface saturation specific humidity is evaluated on a fixed sigma level a = 0.96. The 
free-tropospheric contribution to the saturation specific humidity difference A v q s is evaluated on a sigma level that is a 
fixed amount larger than the sigma level of the tropopausc, a t : a = a t +0.15. The scale estimates P*/A v q s and P*/[q s b] 
are evalu ated within 8° of the equator and (P)/ (q s b) is evaluated using global means. Zonal averages are used for A v q s 
and [q sb ]. 



Equation [TU] assumes that we are averaging over horizontal 
scales that are large enough to include convective down- 
drafts and the associated evaporation of condensate. Ne- 
glecting transients aside from those implicitly included by 
using the truncated vertical velocity, equation (|10[) can be 
decomposed into zonal- means, denoted by [•], and devia- 
tions thereof, denoted by (■)*. The time-mean, denoted by 
(■), of the zonally asymmetric component of the budget is 
then 

- \^d p fo] + W\d P q! + T^d.^ - p^dp^] } sa 

(11) 

Terms 2-4 on the left-hand side, which include zonal asym- 
metries in saturation specific humidity q* , are negligible as 
a result of the weak tem perature gradi e nts in the free tro - 
posphere of the tropics (|Charnevlll963t ISobel et al.ll2001l ). 
In the simulations, these terms combined are an order of 
magnitude smaller than the terms in the dominant balance 



{^*d p [q s }} 



(12) 



Our interest is in the zonally asymmetric component 
of the tropical overturning circulation (denoted w*), the 
Walker circulation, and so we estimate the mass fluxes of 
the ascending or descending branches from (fT"2")) as 



or 
9 



9 



A v q s ' 



(13) 



where A v q s is a zonal- mean vertical saturation specific hu- 
midity difference. 

This is closely related to the thermodynamic budget in 
the tropics where the stratification is close to moist adia- 
batic. On a moist adiabat (at constant saturation equiv- 
alent potential temperature 9*), the vertical gradients of 
potential temperature 9 and saturation specific humidity 
(assuming constant latent heat of vaporization) are related 
by 



(T/9)d p 9\ e « « (L/c p )d p q s \ e * 



(14) 



(e.g., Iribarne and Godson 198ll). Thus, as disc ussed in 
Held and Sodenl (|2006l ) and lSchneider et all (|2010h , the wa- 
ter vapor budget and thermodynamic budget are closely 
linked. The zonally asymmetric component of both bud- 
gets is an approximate two-term balance between vertical 
advection of the zonal-mean potential temperature or satu- 
ration specific humidity by the zonally asymmetric velocity 
an d the zonally asymmetr ic precipitation. 

Held and Soden ( 20061 ) suggested that the free-tropospheric 
saturation specific humidity in scalin g estimates for th e 



upward mass nux is neg ligible or, as m metis ana ma gwai 
(|l989f) and lBettsI (| 19981) . it is linearly related to the bound 
ary layer saturation specific humidity [q s b\- In equation 
(fT3|) . A v q s ~ [q s b] leads to the scaling estimate 



9 



(15) 
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Schneider et al. ( 2010() argued that the free-tropospheric 
contribution is not generally negligible — that is, an actual 
difference in saturation specific humidity must be consid- 
ered in (Tf3| . Physically, this can occur as the result of the 
amplification of warming aloft for moist adiabatic strat- 
ification. However, the amplification of temperature in- 
creases at fixed pressure levels can be offset by increases in 
the pressure depth of the mass fluxes, as occur with a ris- 
ing tropopause. A concrete example of offsetting changes 
in stratification and pressure de pth is given by the "Fixed 
Anvil Temperature" hypothesis ( Hartmann and Larsonl2002 ) . 
which posits that the temperature that convection reaches 
is approximately climate invariant. If this hypothesis holds 
precisely, the free-tropospheric saturation specific humidity 
at the tropical t r opopa use would also be climate invariant. 
Schneider et al. ( 2010l ) also showed that a scaling that in- 
cludes the free-tropospheric contribution to the saturation 
specific humidity difference accounted better for changes in 
the upward component of zonal-mean tropical circulation 
in simulations with zonally symmetric forcing and bound- 
ary conditions. 

Whether or not the free-tropospheric water vapor con- 
tributes depends on the vertical velocity profile Cl{p) that is 
used to separate the integrand: {ui^d p q s } « uj^{£l(p)d p q s }, 
where uft depend s on the ho rizontal, but not the vertical, 



free-tropospheric saturation specific humidity contribution 
to A v q s (blue) and one that neglects it A v q a ~ [q s b] (ma- 
genta). The two are largely indistinguishable and both 
account for the simulated changes in uj*. The scaling that 
neglects the free-tropospheric saturation specific humidity 
is an approximation, so while it appears to be more accu- 
rate in the warm limit, this is the result of compensating 
errors. The scaling constants, obtained by least squares, 
are approximately 1.5 which is consistent with the factor 
of 1.6 that relates uj^* and uj* . 

Figure [5] also shows the global-mean scaling estimate, 
—g(P)/(q s b), in orange. This is similar to the convention 



coordinates (e.g., Sobell 2007 ). This is discussed further in 
the appendix. 

c. Assessment of scaling estimate 

The left panel of Fig. [5] shows that the zonally asym- 
metric component of the pressure velocity in the ascent 
region scales with the zonally asymmetric component of 
the upward pressure velocity: uj* ~ uA* . However, the 
zonally asymmetric total pressure velocity uj* is systemat- 
ically larger (by about 60%) than uj^* . The zonal- mean 
upward vertical velocity [u>^] is larger than the zonal-mean 
total vertical velocity [uj]; hence, when the zonal-mean is 
removed to form the zonally asymmetric component (e.g., 
uj* = uj — [uj]), the total velocity uj* is larger than the up- 
ward velocity a/I"*. The fact that uj* ~ uj^* suggests that 
if the hydrological cycle scaling estimate accounts for the 
changes in the upward pressure velocity, to which it is more 
directly related (|10[) . it will also account for the changes 
in the total pressure velocity, i.e., th e Walker circula t ion. 
However, consistent with the results of lSchneider et al. 
the simulated changes in the zonal-mean pressure velocity 
[uj] in the ascending branch of the Hadley circulation are 
not straightforwardly related to the changes in the upward 
component of the zonal-mean pressure velocity [uj^] . 

The right panel of Fig. [5] shows that the hydrological 
cycle estimate (|13|) applied locally accounts for the changes 
in the Walker circulation. There are two variants of the lo- 
cal scaling in the right panel of Fig. [5] one includes the 



used i n the analysis of IPCC simulat ions in lHeld and Soden 
( 20061 ) and lVecchi and Sodenl (|2007h . who further approxi- 
mated the saturation specific humidity changes by lineariz- 
ing its dependence on temperature about the present cli- 
mate (i.e., 5(q s b) ~ 0.07 x 8(T S }). The global-mean scaling 
captures the gross magnitude of the weakening, but not the 
detail. It tends to underestimate the circulation changes for 
climates colder than the reference climate; at the reference 
climate, the global-mean scaling decreases by 3.1% K _1 , 
compared with uj* which decreases by 4.4% K . 

As the local scaling convention captures the variations 
in the Walker circulation with climate, we can decompose 
it to determine the relative contributions from the precip- 
itation changes (Fig. HI left panel) and the saturation spe- 
cific humidity changes (Fig.|6l right panel). Two factors — 
approximately constant P* and rapidly decreasing A^g^ 1 

or [q s b] — account for the rapid decrease in the Walker 
circulation with warming. 

The zonally asymmetric component of tropical precip- 
itation, P*, is approximately constant with climate; how- 
ever, this is not a general result. Alternative convection 
scheme parameters, in particular relaxing to 90% relative 
humidity instead of 70%, can produce non-monotonic changes 
in P* with climate. The locally averaged version of the 
scaling for Walker circulation changes (fT3|) is similarly ac- 
curate in the alternate set of simulations. These simula- 
tions have increases in P* and a less rapid decrease in the 
Walker circulation uj* with warming from the Earth-like 
reference climate than the simulations presented here. The 
multi-model mean of the models participating in the fourth 
IPCC assessment also features increased P* with warming 



( Meehl et a l. 2007). 



Changes in local precipitation cannot be understood in- 
dependently of the circulation changes as local water vapor 
convergence depends, in part, on the circulation changes. 
In the simulations presented here, the changes in the zonally- 
asymmctric component of net precipitation in the deep 
tropics have a large dynamic component, i.e., a component 
associated with changed winds and fixed specific humidity. 
The dynamic changes reduce the zonally-asymmetric pre- 
cipitation due to weakening surface easterlies with warming 
and offset the thermodynamic component of the changes, 
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Fig. 6. Left: precipitation vs. global-mean surface temperature for zonally asymmetric tropical-mean P* (blue) and 
global-mean (P) (orange). Right: inverse specific humidity vs. global-mean surface temperature for tropical-mean A v q s 
(blue), tropical-mean [q s i,} (magenta), and global- mean {q s b) (orange). The averaging conventions are the same as in 
Fig. EJ 



i.e., the component associated with changes in specific hu- 
midity and fixed winds, which increase the the zonally- 
asymmetric precipitation with warming. This is qualita- 
tively consistent with the multi-model mean of th e fourth 



that consists of a two-term balance between the zonally 
asymmetric component of the Q-flux and evaporation. It 
leads to a scaling for the zonal difference in tropical surface 
temp erature that depends inversely on the rate of change 



IPCC assessment simulations presented in Fig. 7 of lHeld and 
( 20061 ): the equatorial Pacific is a region in which the 
changes in net precipitation arc quite different from a ther- 
modynamic estimate that assumes fixed circulation. 

The saturation specific humidity is a rapidly increasing 
function of temperature, so, as expected, its inverse rapidly 
decreases with warming. The difference between the in- 
verse saturation specific humidity when the free-tropospheric 
component of A v q s is neglected and when it is retained is 
small. This is a consequence of the depth of the circulation — 
it reaches sufficiently low pressures and temperatures that 
q s is negligible in that region of the free troposphere. The 
appendix has a more extensive discussion of the depen- 
dence of the scaling estimate on the vertical structure of 
the circulation and how the results of the simulations pre- 
sented here relate to previous work. 

Figure [6] also sheds light on why the global-mean con- 
vention may, in some cases, appear adequate. Global-mean 
precipitation (P) increases more rapidly than the zonally 
asymmetric component of tropical precipitation P* ; global- 
mean surface saturation specific humidity (q s b) increases 
more rapidly than the tropical-mean surface saturation spe- 
cific humidity difference A„g s . As the scaling depends on 
the ratio, the more rapid changes in the individual com- 
ponents largely cancel, but this must be regarded as fortu- 
itous. 

5. Conclusions 

Zonal surface temperature gradients decrease rapidly as 
the climate warms in the idealized GCM. This can be un- 
derstood by considering a simplified surface energy budget 



::mp 

Soafera turation specific humidity (or the Clausius-Clapeyron 



relation) , evaluated using the mean surface temperature in 
the tropics if the Q-flux is fixed. 

The Walker circulation rapidly weakens with warming 
in the idealized GCM. The changes can be accounted for 
using locally evaluated hydrological cycle scaling estimates. 
The changes in the scalings can be decomposed into i) equa- 
torial precipitation increases being approximately zonally 
uniform with warming and ii) rapid, Clausius-Clapeyron 
increases in near-surface saturation specific humidity with 
warming. The global-mean version of the hydrological cy- 
cle scale estimate does not accurately account for the changes 
in the Walker circulation in all climates, though it does 
capture gross changes due to compensating changes in pre- 
cipitation and saturation specific humidity that are larger 
than the locally averaged quantities. 

That the hydrological cycle scaling estimate requires 
knowledge of local precipitation changes is a limitation. 
Ideally, one would want to understand circulation changes 
independently and use that result in the water vapor bud- 
get to constrain local precipitation changes. The moist 
static energy bud get is one ap proach, though it has estab- 
lished difficulties ISobell 120071 ) . Another approach is tha t 
of the Lindzen-Nigam model ([Lindzen and Nigaml II9871 ). 
in which the surface temperature field is used to deter- 
mine pressure gradients, which in combination with a lin- 
ear momentum equation determine boundary layer conver- 
gence. We have solved the Lindzen-Nigam equations us- 
ing the GCM simulated surface temperature fields as forc- 
ing and found that the changes the Lindzen-Nigam model 
predicted in the boundary layer mass flux underestimated 
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those simulated by the GCM by roughly a factor of two over 
the range of climates. Devising conceptual models of how 
tropical circulations change with climate that adequately 
account for GCM simulations remains a challenge. 

A limitation of this study is the prescribed and invariant 
Q-flux. But understanding the atmosphere-only dynamics 
is a necessary step towards to the more realistic case that 
includes an interactive ocean. 
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APPENDIX 



Dependence of water vapor scaling on vertical 
structure 

To understand why the frec-troposphcric contribution 
to the saturation specific humidity difference is negligi- 
ble, we examine the vertical structure of the vertical ve- 
locity, as A v q s is an approximation of {tt(p)d p q s }. Also, 
the vertical structure of the circulation has importance be- 
yond the hydrological cycle scaling estimate: it determines 
whether moist static en ergy is imported or exported by 
overturning circulations (jSobel 2007 ; Back and Bretherton 
l2006t IPeters et al] l2008h . If the vertical velocity profile 
has larger amplitude in the free troposphere, it makes the 
saturation specific humidity difference less sensitive to the 
free-tropospheric values entering A v q s because tempera- 
ture and saturation specific humidities arc small there. 
Indeed, the vertical velocity profile is "top-heavy" in the 
ascending branch of the Walker circulation (solid lines in 
Fig. [7|); that is, it has significant amplitude in the upper 
troposphere. In contrast, the zonal-mean vertical velocity 
profile has more amplitude near the surface and less ampli- 
tude in the upper troposphere (dashed lines in Fig. [TJ). The 
vertical velocity profiles of the idealized simulations are 
broadly consistent with the regional variations in the verti- 
cal ve locity profile in Earth's tropics (jBack and Bretherton 
l2006h . 



For the zonal-mean component of the simulations pre- 
sented here, the hydrological cycle scaling estimate that 
neglects the changes in the saturation specific humidity 
away from the boundary layer (A v q s ~ [q s b]) is distinguish- 
able from the scaling estimate that retains the full A v q s 
(i.e., the representative pressure levels when the integral 
is discretized are closer to each other than in the zonally 
asymmetric case and the free-tropospheric contribution is 
not negli gible). This is consis tent with the results pre- 



sented in ISchneider et ah I <l2010h 



Figure [7] also shows that the pressure depth over which 
there is significant vertical velocity amplitude increases and 
the tropopause pressure decreases with warming. For this 
reason, we vary the pressure at which the free-tropospheric 
contribution to A v q s is evaluated in Fig.[SJ neglecting these 
changes can potentially lead to different scaling behavior. 

There are two ways in which the vertical velocity profile 
can lead to Clausius-Clapcyron-like scaling of A v q s . First, 
a larger pressure depth decreases the free troposphere's 
weight (e.g., in the limit of the whole troposphere, the 
tropopause saturation specific humidity is orders of mag- 
nitude smaller than that of the surface). Second, hold- 
ing the pressure depth fixed and shifting it upward (e.g., 
the difference between 800 hPa and 300 hPa compared to 
the difference between 900 hPa and 400 hPa) leads to more 
rapid changes in the saturation specific humidity difference 
as the temperatures are lower and the rate of change of the 
saturation specific humidity is larger. These changes will 
appear closer to the rate given by the Clausius-Clapcyron 
relation if it is evaluated at a fixed level such as the surface 
(related issues about the sensitivity of saturation specific 
humidity to where it is evaluated and the possible impre- 
cision that results are discussed in lO'Gorman and Muller 
20101) . 



In an analysis of comprehensive GCMs. lVecchi and Soden 



(|2007j ) showed that changes in the global-mean upward 
component of the vertical velocity scaled with the hydrolog- 
ical cycle scaling estimate that neglects the free-tropospheric 
contribution to the saturation specific humidity (i.e., (ur) ~ —g (P) / (i 
but the scaling overestimated the magnitude of the vertical 
velocity ch anges by about a f a ctor o f two to three. How- 
ever, when IVecchi and Sodenl ( 20071 ) considered only up- 



ward pressure velocities exceeding 0.05 Pas -1 (i.e., ur = 
ojH(— 0.05 Pas -1 — ui)), the changes in the vertical veloc- 
ity were larger and consistent with the scaling estimate. 
Changes in the vertical structure of the circulation may ex - 
plain why the scaling in Fig. 5 of IVecchi and Sodenl (|2007l ) 
improved when a non-zero threshold for upward vertical 
velocity was introduced: With a threshold of zero, the ba- 
sis function is less top-heavy and the changes in the free- 
tropospheric saturation specific humidity contribution to 
A v q s cannot be neglected, in which case approximating 
A v q s with [q s b] leads to an overestimate of the predicted 
weakening of the circulation with warming; however, for re- 
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Fig. 7. Vertical structure of vertical velocity f2(cr) = 
w*/min(o;^) for simulations with global-mean surface tem- 
perature (T s ) 289.9 K (black) and 295.3 K (gray). Solid 
lines are zonally asymmetric pressure velocity profiles fi* 
evaluated in the ascending region (cf. Fig. [T|) and dashed 
lines are zonal-mean pressure velocity profiles [f2]. The 
dashed-dotted line is the zonal-mean tropical tropopause 
defined by the WMO 2Kkm~ 1 lapse rate criterion. 

gions of strong ascent (such as those considered here), the 
vertical velocity profile can be sufficiently deep that the 
free-tropospheric contribution to A v q s is negligible, A v q s 
scales similarly to [q s b] , and the circulati on weakens rapidly 
with w arming. Another factor is that Vecchi and Sodenl 
( 20071 ) changed from a threshold of Pas -1 , which includes 
extratropical regions of ascending motion, to a threshold 
of 0.05 Pas -1 , which limits the area considered to a few, 
strongly ascending tropical locations, while using global- 
mean precipitation and surface specific humidity in both 
cases. The utility of global-mean scaling estimates for lo- 
cal circulation changes is assessed in Fig.[5]and Fig.|5]which 
show that it does not generally account for the simulation 
results. Because changing the threshold changes both the 
geographic area under consideration and vertical structure, 
it is unclear which is more important. 
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